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Abstract 

In an adiabatic rapid passage experiment, the Bloch vector of a two-level system (qubit) is 
inverted by slowly inverting an external field to which it is coupled, and along which it is initially 
aligned. In twisted rapid passage, the external field is allowed to twist around its initial direction 
with azimuthal angle (j){t) at the same time that it is inverted. For polynomial twist: (f){t) ~ 
Bt^. We show that for n > 3, multiple avoided crossings can occur during the inversion of the 
external field, and that these crossings give rise to strong interference effects in the qubit transition 
probability. The transition probability is found to be a function of the twist strength B, which 
can be used to control the time-separation of the avoided crossings, and hence the character of the 
interference. Constructive and destructive interference are possible. The interference effects are a 
consequence of the temporal phase coherence of the wavefunction. The ability to vary this coherence 
by varying the temporal separation of the avoided crossings renders twisted rapid passage with 
adjustable twist strength into a temporal interferometer through which qubit transitions can be 
greatly enhanced or suppressed. Possible application of this interference mechanism to construction 
of fast fault-tolerant quantum CNOT and NOT gates is discussed. 

PACS numbers: 03.67.Lx,07.60.Ly,31.50.Gh 
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I. INTRODUCTION 



Adiabatic rapid passage fARP) is a well-known procedure for inverting the Bloch vector of 
a two-level system (qubit) This is accomplished by inverting an external field F{t) which 
couples to the qubit, and along which the qubit is initially aligned. The field inversion is done 
on a time-scale that is large compared to the inverse Rabi frequency uJq'^ (viz. adiabatic), 
though small compared to the thermal relaxation time r (viz. rapid). In the usual case, F{t) 
remains within a plane that includes the origin: F(t) = 6x + at z, with —To/2 <t< To/2, 
and uJq^ ^ Tq ^ r. 

ARP can be used to implement a NOT gate on the quantum state of a qubit. If one 
identifies the computational basis states |0) and respectively, with the spin-up and spin- 
down eigenstates along the initial direction of the external field F(— To/2), then ARP maps 
|0) ^ |1) which is the defining operation of a NOT gate. Occurrence of a transition during 
ARP corresponds to an error in the NOT gate since the Bloch vector is not inverted, and 
thus \i) l^), {i = 0,1). Thus we can identify the ARP transition probability with the 
NOT gate error probability. For ARP, the adiabatic nature of the inversion ensures that 
the error probability is exponentially small. The price paid for this reliability, however, is 
an extremely slow NOT gate. 

In twisted adiabatic rapid passage, the external field is allowed to twist around its ini- 
tial direction with azimuthal angle 0(t) at the same time that it is adiabatically inverted: 
F{t) = b cos (j){t) X + 6 sin (f){t) y + atz. Reference showed that for twisted ARP, the expo- 
nentially small transition probability contains a factor exp[rg] of purely geometric origin. 
The simplest case where 7^ corresponds to quadratic twist: 0(t) = Bt^. Zwanziger et. 
al. 01 were able to experimentally realize ARP with quadratic twist and obtained results in 
agreement with the predictions of Reference Q. 

In this paper we will consider twisted rapid passage with polynomial twist, (j){t) ~ Bt^, 
and we will focus exclusively on qubit inversions done at non-adiahatic rates. Although we 
will briefly consider quadratic twist in Section |n] as a test case for our numerical simulations, 
our interest will not be the geometric effect of Reference Q. Instead, our primary focus will 
be on establishing the existence of multiple avoided crossings during twisted rapid passage 
when n > 3, and with exploring some of their consequences. After general considerations 
(Section |n|), we will explicitly examine cubic (n = 3) and quartic (n = 4) twist, and will 



provide clear evidence that the multiple avoided crossings produce strong interference effects 
in the qubit transition probability. The transition probability is shown to be a function of the 
twist strength B, which can be used to control the time-separation of the avoided crossings, 
and hence the character of the interference (constructive or destructive). Cubic and quartic 
twist are examined in Sections IIIII and IIV| respectively. We shall see that interference 
between the multiple avoided crossings can greatly enhance or suppress qubit transitions. 
The interference effects are a direct consequence of the temporal phase coherence of the 
wavefunction. The ability to vary this coherence by varying the temporal separation of 
the avoided crossings renders twisted rapid passage with adjustable twist strength into a 
temporal interferometer through which qubit transitions can be controlled. It will be shown 
that quartic twist can implement qubit inversion non-adiabatically while operating at a 
fidelity that exceeds the threshold for fault tolerant operation. Finally, in Section |3 we 
summarize our results and discuss possible application of this interference mechanism to 
quantum computing. In particular, we describe how one might use non-adiabatic rapid 
passage with quartic twist to construct a fast fault-tolerant quantum CNOT gate. 

It is worth noting that experimental confirmation of the work described in this paper has 
recently been carried out by Zwanziger et. al. [4]. They have realized non-adiabatic rapid 
passage with both cubic and quartic twist, and have observed clear evidence of constructive 
and destructive interference in the qubit transition probability due to interference between 
the avoided crossings, with excellent agreement between the experimental data and our 
numerical simulations. This experimental work provides clear proof-of-principle for our 
thesis that controllable quantum interference exists during twisted rapid passage. With this 
thesis now experimentally confirmed, future research can focus on applying this interference 
to the task of constructing fast fault tolerant quantum CNOT and NOT gates. 

After this paper was submitted, previous work was brought to our attention which also 
examined models of rapid passage in which more than one avoided crossing is possible, and 
in which interference effects were also considered ^ 0, llj. These papers focused solely on 
the adiabatic limit. Application of this adiabatic theory to the Zwanziger experiment yields 
predictions that are in poor agreement with the experimental results [4]. This failing is no 
doubt a consequence of the non-adiabatic character of this experiment whose results thus 
lies beyond the scope of the adiabatic theory developed in these papers. In contrast, the 
work we present below is principally interested in the non-adiabatic limit, and our simulation 



results are in full agreement with experiment We also consider possible application of 
these interference effects to the construction of fast fault-tolerant quantum CNOT and NOT 
gates. Ref. do not consider such applications. 

II. TWISTED RAPID PASSAGE 

We begin by briefly summarizing the essential features of rapid passage in the absence 
of twist. Twistless rapid passage describes a wide variety of phenomena, ranging from 
magnetization reversal in NMR, to electronic transition during a slow atomic collision. The 
essential situation is that of a qubit which is Zeeman-coupled to a background field F(t), 
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H{t) = a- . Fit) = , (1) 

\b -at I 

with F{t) = 6x + atz. This particular form for F(t) describes inversion of the background 
field in such a way that it remains in the x-z plane throughout the inversion. For simplicity, 
we assume a, b > throughout this paper. The instantaneous energies E±{t) are: 



E^{t) = ±^WTW , (2) 

and an avoided crossing is seen to occur at t = where the energy gap is minimum. The 
Schrodingei_cl^namics for twistless rapid passage can be solved exactly for arbitrary values 
, and yields the Landau- Zener expression for the transition probability Plz'- 



of a and b M 



Plz = exp 



TTb' 



(3) 



h\a\ 

A. Twisted Rapid Passage and Multiple Avoided Crossings 

In twisted rapid passage, the background field F{t) is allowed to twist around its initial 
direction during the course of its inversion: F{t) = b cos x + 6 sin y + at z. It proves 
convenient to transform to the rotating frame in which the x-y component of the background 
field is instantaneously at rest. This is accomplished via the unitary transformation U{t) = 
exp[—{i/2)(j){t)(Jz]. The Hamiltonian H{t) in this frame is: 

_ at - |0 ) b 



where F{t) = 6x + (at — ^0/2) z is the background field as seen in the rotating frame, 
and a dot over a symbol represents the time derivative of that symbol. The instantaneous 
energy eigenvalues are E±{t) = ±^ (^at — (/i0)/2^ + 6^. Avoided crossings occur when 
the energy gap is minimum, corresponding to when 

h dd) „ 

at r = . (5) 

2 dt ^ ^ 

For polynomial twist: </)„(t) = c„i?t", where B is the twist strength. The dimensionless con- 
stant c„ has been introduced to simplify some of the formulas below. For later convenience, 
we chose c„ = 2/n. For polynomial twist, it is easily checked that eqn. ^ always has the 
root: 

t = , (6) 
and that for n > 3, eqn. (jSj) also has the n — 2 roots: 



All together, equation ^ has n — 1 roots, though only the real roots correspond to avoided 
crossings. For quadratic twist {n = 2), only eqn. (jH)) arises. Thus, for this case, only the 
avoided crossing at t = is possible. For n > 3, along with the avoided crossing at t = 0, real 
solutions to eqn. (|7j) also occur. The different possibilities for this situation are summarized 
in Table m We see that for polynomial twist with n > 3, multiple avoided crossings always 

TABLE I: Classification of regimes under which multiple avoided crossings occur for polynomial 
twist with n > 3. 

1. sgnB = +1 

(a) n odd; 2 avoided crossings at: t = and t = (a/hB)^^^^"'^^ 

(b) n even; 3 avoided crossings at: t = and t = it {a/hB)^^^^~'^^ 

2. sgnB = -1 

(a) n odd; 2 avoided crossings at: t = and t = — {a/h\B\)^^^^^'^^ 

(b) n even; 1 avoided crossing at: t = 



occur for positive twist strength B, while for negative twist strength, multiple avoided 
crossings only occur when n is odd. Note that the time separating the multiple avoided 
crossings can be adjusted by variation of the twist strength B and/or the inversion rate a. 



B. Brief Detour: Background on Quadratic Twist 

As mentioned earlier, quadratic twist has already been examined in the literature P|. It 
is of interest here only because its dynamics can be solved exactly, and thus allows us to test 
our numerical simulations before proceeding to unexplored cases of twisted rapid passage. 
For quadratic twist, 02 = 2Bt. Inserting this into eqn. ^ gives F(t) = b'k + atz, with 
a = a — h\B\ (sgni?). Thus rapid passage with quadratic twist maps onto twistless rapid 
passage with a ^ a. This allows us to obtain an exact result for the transition probability 
P2 for arbitrary values of a and b from eqn. (jS)) with a — > a: 



exp 



(8) 



h\ a — h\B\ (sgni?) 

In the adiabatic limit, this reduces to P2 = Plz ^w[^g]j where Tg = —iiBb'^la? is the 
geometric exponent discovered in Ref. Eqn. (jS)) makes the interesting prediction that 
a complete quenching of transitions will occur when sgni? = +1 and a = hB, while no 
such quenching is possible for sgni? = — 1. Zwanziger et. al. |3| were able to realize rapid 
passage with quadratic twist experimentally and confirmed the existence of Tg, and the 
twist-dependent quenching of transitions. We now show that our numerical simulation also 
reproduces these effects. 

C. Simulation Details 

The equations that drive the numerical simulation follow from the Schrodinger equation 
in the non-rotating frame: 

^n^^\^P)=H{t)\^) , (9) 

where H{t) = cr-F(t), and F(t) = 6cos0(t) x + 6sin</)(t) y + at z. To obtain these equations 
in the adiabatic representation, we expand I'lpit)) in the instantaneous eigenstates \E±(t)) 
ofH{ty. 



Here 7±(t) are the geometric phases [10| associated with the energy- levels E±{t), respectively, 
and 

^^{t) = z{E^{t)\^^\E^{t))=z{E^{t)\E^{t)) . (11) 



Substituting eqn. (fUH) into 0, and using the orthonormality of the instantaneous eigen- 
states, one obtains the equations of motion for the expansion coefficients S{t) and I{t): 



^ = -r(t) e-'-^-^o/2^^^(^) /(t) , (12a) 
^ = r(t) e^Z-To/^^ieW s{t) . (12b) 

(J/V 



Here, 



m = ^i^^^^^-(74t)-7-(t)) , (13) 
m = {E4t)\E4t)) , (14) 

and one can show that r*{t) = — {E^{t) \ E^{t) ). Eqns. (fT^ are the qubit equations of mo- 
tion in the adiabatic representation and include the influence of the geometric phase on the 
dynamics through S{t). In the case of twistless rapid passage, the geometric phase vanishes, 
and eqns. ()12|) reduce to the well-known equations of motion for a two-level system found 
in Ref. Eqns. (|T2|l can be put in dimensionless form if we introduce the dimensionless 
variables: r = {a/b)t, T = {b/a)r, and 6 = {b/a)6. Here a and b are the parameters that 
appear in the background field F{t). One obtains: 

= -r e-^^-^o/^^^^W/(r) , (15a) 

ar 

^ = Te'^-^o/^'''^'^ Sir) , (15b) 

where tq = (a/6)To is the (dimensionless) time over which the qubit evolves. For rapid 
passage, the qubit is initially in the negative energy level \E_{—tq/2) ). This corresponds to 
the initial condition: 

S{-To/2) = 1 , (16a) 
/(-ro/2) = . (16b) 

Our numerical simulation integrates eqns. (fTHjl over the time- interval [—ro/2, ro/2] subject 
to initial condition ()16|). From this we determine the asymptotic transition probability P: 

P=\I{ro/2)\' , (17) 



for To ^ 1. Later, we will need the r- values corresponding to the avoided crossings. These 
are determined by rewriting eqns. (jH)) and ((Tj) in dimensionless form. To this end, we 
introduce 

Vn = —-r- , (18) 



and recalling that r = {a/b)t, one easily obtains: 



r = 



(19) 



and 



r = [sgnr], 



(20) 



The avoided crossings correspond to r = and also, for n > 3, the real solutions of eqn. ((201) ■ 



D. Simulation Test Case: Quadratic Twist 

For quadratic twist 02 (t) = Bt"^. The instantaneous eigenvalues and eigenvectors of H{t) 
are easily found to be E±{t) = ±E{t), where E{t) = ^JW' + (at)^, and 

cos % \ I sin^ , 

\E+it))= ' ; \E-{t))= ' (21) 

sin|e''^2 ; I _cos|e''^2 
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with cos 6' = at/E. From the eigenstates one obtains: 

i±{t) = -^(iTcos^) ; (22a) 

m = l-^^sine ; (22b) 
2E 

6(t) = — -02cos^ . (22c) 
n 

r(r) and 6{t) are then determined from eqns. (j22b() and ()22cj) and are found to depend 
parametrically on the dimensionless "inversion rate" A = ha/b'^ and the dimensionless "twist 
strength" ri2 = hB/a. "Inversion rate" and "twist strength" are placed in quotes as A 
does not depend solely on the inversion rate a, nor ri2 solely on the twist strength B. 
Crudely speaking, A = 1 can be thought of as the boundary separating adiabatic and non- 
adiabatic inversion rates, with A > 1 corresponding to non-adiabatic inversion. Having 
determined r(r) and 5(r), eqns. (fT^ are integrated numerically using an adjustable step- 
size fourth-order Runge-Kutta algorithm. To simplify comparison of the numerical result 



for the transition probability with the exact result P21 we re- write eqn. 
7] 2- One finds: 

P2 = exp 



vr 



1 



A \ l-r]2 

Figure H shows a representative plot of the transition probability P(t) 

Transition Probability vs. Time 



in terms of A and 



(23) 

|/(r)P versus r. It 




FIG. 1: Representative plot of transition probability P{t) for quadratic twist with A = 10.0 and 
r?2 = 2.5. 



is clear for the Figure that the transition occurs in the vicinity of the avoided crossing at r = 
0. Note also that P{t) has a small oscillation about its asymptotic value P = lim.r-+oo P{t)- 
To average out the oscillation, P{t) (for given A and 772) was calculated for 10 different 
values of r ^ 1, and P was identified with the average. Figures |21 and 121 show our numerical 
results for P for various values of for A = 10.0 and A = 3.0, respectively. Also plotted 
in each of these Figures is the exact result P2 (eqn. ()23|) ). Figures |2 and El show that our 
numerical results are in excellent agreement with the exact result P2; and clearly show the 
quenching of transitions at 772 = 1, and the absence of quenching for negative ri2- The A 
values shown are purposely highly non-adiabatic. We see that the twist-induced quenching 
clearly persists into the non-adiabatic regime, although the width of the quench decreases 
with increasing A. The agreement of our simulations with eqn. (|^^ at small 772 indicates 
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FIG. 2: Numerical results for the asymptotic transition probability P versus r]2 for quadratic twist 
with A = 10.0. Also plotted is the exact result P2- 

that our simulations also account for the geometric factor exp [Tg] in P2. Having established 
that our numerical algorithm correctly reproduces the essential results of rapid passage with 
quadratic twist, we go on to consider the unexplored areas of rapid passage with higher order 
twist. Referring to Table HI we see that all cases with odd n have 2 avoided crossings. Cubic 
(n = 3) twist corresponds to the simplest example of odd-order twist, and it is examined in 
the following Section. Similarly, quartic (n = 4) twist is the simplest example of even-order 
twist, and we examine it in Section HVl 

III. CUBIC TWIST 

Having successfully tested our numerical algorithm against the exact results for quadratic 
twist, we go on to consider cubic twist for which 03(t) = {2/3)Bt^, and rj^ = KBh/a^ (see 
eqn. p8j) ). As in Section |H1 the instantaneous eigenvalues of H{t) are E±{t) = ±E{t), 
and the instantaneous eigenstates are given by eqn. (PT|) with 02 (i) (psit)- Eqns. 
again apply, however 02 (^) 'Psit), and r(r) and 6{t) are determined from eqns. ()22b|) and 
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FIG. 3: Numerical results for the asymptotic transition probability P versus r]2 for quadratic twist 
with A = 3.0. Also plotted is the exact result P2- 

fl22c|) . Having determined r(r) and 6{t), eqns. ()15|) can be numerically integrated subject 
to the initial condition specified in eqns. ()lfi|l . Before examining results of that integration, 
we show in Figure 0] a plot of the numerical results for the transition probability P(r) for 
A = 5.0 and r^s = 0. This corresponds to twistless non-adiabatic rapid passage, and we 
include this plot for later comparison with related plots for cubic and quartic twist. The 
asymptotic transition probability for this case is P = 0.533. Thus, if we were to use this 
example of twistless non-adiabatic rapid passage to implement a fast NOT-operation on a 
qubit, the operation would be slightly more likely to produce an inversion (bit-flip) error 
than not. We will show below that if a small amount of cubic twist is included, the bit-flip 
error probability can be reduced by 2 orders of magnitude while still maintaining the non- 
adiabatic inversion rate A = 5.0. This substantial reduction in error probability is due to 
destructive interference between the two avoided crossings that occur during rapid passage 
with cubic twist. 
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FIG. 4: Plot of the transition probability P{t) for twistless non-adiabatic rapid passage with 
A = 5.0 and 773 = 0. 

A. Demonstration of Quantum Interference 

From eqns. (fTT^ and (pUjl . we see that cubic twist is expected to have 2 avoided crossings 
at Ti = and T2 = sgnrj^/lrisl. Figures El and El show P(r) for A = 5.0 and 773 = 0.02 
and ?73 = —0.02, respectively. Figure El® clearly shows the expected avoided crossings at 
r = and r = 50 (—50). It is also clear from these Figures, and comparison with Figure |5 
that the avoided crossings are constructively interfering, leading to an asymptotic transition 
probability of P = 0.997. Figures O and El show P{t) for A = 5.0 and 173 = 0.05 and —0.05, 
respectively. The avoided crossings in Figure (jHl) clearly occur at r = and r = 20 
(—20) as expected. Here the avoided crossings interfere destructively, with P = 0.270. 
Summarizing, we see that: (1) two avoided crossings do occur during rapid passage with 
cubic twist as predicted in Table U (2) the avoided crossings produce interference effects in 
the asymptotic transition probability P which can be controlled through variation of their 
separation; and (3) the separation of the avoided crossings Ar^c = |t2 — ti| = 1/|773| can 
be altered by varying rj^ = hBh/a?. We now consider two possible applications of this 
interference effect. 
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FIG. 5: The transition probability P{t) for non-adiabatic rapid passage with cubic twist with 
A = 5.0 and 773 = 0.02. 

B. Non-Resonant Pump 



First, consider twistless adiabatic rapid passage with A = 0.5 and 773 = 0. Figure IHl 
show the transition probabihty P{t) for this case. The asymptotic transition probabihty is 
P = 1.87 X 10"'^. Figure ITUl shows P{t) for adiabatic rapid passage with cubic twist with 
A = 0.5 and rj^ = 0.04. The asymptotic transition probabihty in this case is P = 0.996! 
Thus, by introducing a small amount of cubic twist, constructive interference between the 
avoided crossings transforms adiabatic rapid passage into a non-resonant pump for the qubit 
energy levels. Figures El and El indicate that, should it be desired, equally large transition 
probabilities are also possible at faster inversion rates A. It is worth pointing out that to 
produce such a large transition probability using twistless non-adiabatic rapid passage would 
require A = 784 (see eqn. with 772 = 0) as opposed to A ~ 0.5 — 5.0 when cubic twist is 
exploited. 
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FIG. 6: The transition probability P{t) for non-adiabatic rapid passage with cubic twist with 
A = 5.0 and 773 = -0.02. 

C. Transition Quenching 

We now show that one can utihze the interference between avoided crossings to strongly 
suppress qubit transitions during non-adiabatic rapid passage with cubic twist. Figure ITTl 
shows -P(t) for A = 5.0 and rj^ = 4.577 x 10^^. The asymptotic transition probability for this 
case is P = 3.44 x 10^'^. This is to be compared with twistless rapid passage with A = 5.0 
(Figure 13)) for which P = 0.533. Destructive interference between the two avoided crossings 
has reduced the transition probability P by 2 orders of magnitude relative to the twistless 
case shown in Figure El Thus if we were to implement a fast NOT-operation using non- 
adiabatic rapid passage with cubic twist at A = 5.0 and r^s = 4.577 x 10"^, we would obtain 
(on average) 1 bit-flip error per 291 NOT-operations. By comparison, twistless rapid passage 
with A = 5.0 would produce (on average) 1 bit-flip error for every 2 NOT-operations. This 
result strongly suggest the value of exploring whether this destructive interference between 
avoided crossings during twisted rapid passage could be exploited to produce fast rehable 
quantum NOT and CNOT logic gates. As striking as this result for cubic twist is, we shall 
see in the following Section that quartic twist can reduce the bit-flip error probability even 
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FIG. 7: The transition probability P{t) for non-adiabatic rapid passage with cubic twist with 
A = 5.0 and 773 = 0.05. 

more dramatically. 



IV. QUARTIC TWIST 

For quartic twist 04(t) = (l/2)i?t^ and rj^ = hBb^/a^. Avoided crossings are expected 
to occur at ti = 0, and at T2 = ±l/y^ (when sgnr]^ = +1; see eqn. and Table HI). 
Formally, the analysis of quartic twist parallels that of quadratic and cubic twist. With the 
substitution 02 (i) (piit), eqns. (PT|) and continue to apply, and one determines r(r) 



and 6{t) from eqns. ()22b|) and ()22c|) . Once r(r) and 6{t) are known, eqns. (|T3j) can be 
integrated numerically subject to the initial condition specified in eqns. f|TH| . 



A. Demonstration of Quantum Interference 

In Figure IT^ we plot the transition probability P{t) for A = 5.0 and rj^ = 4.6 x 10~^. The 
expected avoided crossings at ri = and T2 = ±46.63 are clearly visible. The asymptotic 
transition probability for this case is P = 0.88. For twistless rapid passage with A = 5.0 (see 
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FIG. 8: The transition probability P{t) for non-adiabatic rapid passage with cubic twist with 
A = 5.0 and 773 = -0.05. 

Figure Ej), P = 0.533. Thus the avoided crossings in Figure IT^ are constructively interfering, 
leading to an enhancement of the transition probability P. Figure IT!?1 shows P(t) for quartic 
twist with A = 5.0 and 774 = —4.6 x 10~^. This Figure clearly shows only one avoided 
crossing at ti = 0, as expected for sgnrj^ = —1 (see Table H}. The asymptotic transition 
probability in this case is P = 0.533 which equals the result for twistless rapid passage with 
A = 5.0 (Figure IH) to the level of precision obtained in our calculation. 

Figure IT^ plots P{t) for A = 5.0 and 174 = 1.6 x 10~^. The Figure clearly shows the 
expected crossings at ri = and T2 = ±25.0. The asymptotic transition probability is 
P = 6.93 X 10~^ and corresponds to destructive interference relative to twistless rapid passage 
with A = 5.0 (FigurelH). We do not include a plot of P{t) for A = 5.0 and rj^ = —1.6 x 10~^ 
as it is similar to Figure IT^ one avoided crossing at ri = and P = 0.533. 

Summarizing these results, we see that: (i) three (one) avoided crossings (crossing) oc- 
cur(s) as predicted in Table HI when sgnrj^ = +1 (—1); (ii) the avoided crossings produce 
interference effects in the transition probability, with the character of the interference (con- 
structive or destructive) determined by the separation of the avoided crossings; and (iii) the 
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FIG. 9: The transition probability P{t) for twistless adiabatic rapid passage with A = 0.5 and 
r/s = 0. Note the greatly reduced vertical scale compared to previous figures. 



separation of adjacent avoided crossings is given by Ar^c = |t2 — ti| = 1/ ^/rfl {sgnrj^ = +1), 
and it is controllable through variation of 774 = hBb^/a^. 



B. Non-Resonant Pump 

Quartic twist does not appear to be as effective at pumping the qubit energy-levels as 
cubic twist. Figure Uni shows P{t) for A = 0.5 and 774 = 6.45 x 10^^. The expected avoided 
crossings at ri = and T2 = ±12.45 are clearly visible, and the asymptotic transition 
probability is P = 0.20. Although this is a 2 order of magnitude improvement over twistless 
adiabatic rapid passage with A = 0.5 (FigurelHI), it falls well short of the transition probability 
P = 0.996 easily obtainable with cubic twist. In fact, for 774 < 1, P ~ 0.20 was among the 
largest P-values we could find. If larger values of twist strength are allowed, the largest 
transition probability we could find was P = 0.64 at rj^ = 3.00. 
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FIG. 10: The transition probability P{t) for adiabatic rapid passage with cubic twist with A = 0.5 
and r]3 = 0.04. 

C. Transition Quenching 



Quartic twist proves to be much more effective at quenching transitions during non- 
adiabatic rapid passage than cubic twist. Table E] gives the transition probabilities for 
quartic twist pulses for which A = 5.00 and 174 lies in the interval [3.95 x 10"'^, 4.04 x 10"'^]. 
The essential thing to notice in Table ITTl is that for 774 = 4.00 x 10^^, the transition probability 
P = 4 X 10~^. This is significant for the following reason. It has been shown that a quantum 
computation of arbitrarily long duration becomes possible if the quantum logic gates used 
to implement the computation all have error probabilities (per gate operation) which lie 
below the threshold Pft for fault tolerant operation . This threshold has been estimated 
to be Pft ~ 10^'^ — 10~^ In terms of the gate fidelity F = 1 — P, the more optimistic 
estimate for Pft gives Fft = 0.9999. We see that for A = 5.00 and 774 = 4.00 x 10"^, twisted 
rapid passage with quartic twist gives a gate fidelity of F = 0.99996 which exceeds the best 
case estimate for fault tolerant operation Fft = 0.9999. This fault tolerant performance is 
achieved while inverting the qubit at a non- adiabatic rate. The reader should note that the 
values A = 5.00 and 774 = 4.0 x 10^'^ can be realized with existing NMR technology (see 
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FIG. 11: The transition probability P{t) for non-adiabatic rapid passage with cubic twist with 
A = 5.0 and 7/3 = 4.577 x 10~^. Note the slightly reduced vertical scale. 



Section IV D|) . Our analysis raises the exciting possibility that non-adiabatic rapid passage 
with quartic twist might provide a means of realizing fast fault-tolerant NOT and CNOT 
gates. The novelty of this prospect is the marriage of operational speed with fault-tolerance. 
This marriage of speed and reliability is a direct consequence of the destructive interference 
which is possible between the 3 avoided crossings that arise during rapid passage with 
quartic twist. Quantum CNOT gates are ubiquitous in quantum computing and quantum 
error correction Q|. Thus, determining how to implement them in a fast fault 

tolerant manner is a potentially significant development for the field. 



V. DISCUSSION 
A. Summary 

It has been our aim in this paper to show that multiple avoided crossings can arise 
during twisted rapid passage, and that by varying their time-separation, interference effects 
are produced which allow for a direct control over qubit transitions. This time-separation is 
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FIG. 12: The transition probability P{t) for non-adiabatic rapid passage with quartic twist with 
A = 5.0 and 774 = 4.6 x 10"^. 

controlled through the (dimensionless) twist strength 77, and the resulting interference can 
be constructive (enhancing transitions) or destructive (reducing transitions). For nth-order 
polynomial twist, //„ = hBb'"'~'^ / a^~^ , where B is the (dimensionful) twist strength, 2b is the 
energy-gap separating the qubit energy-levels at an avoided crossing, and a is the inversion 
rate of the external field F(t) (see Section |n}. The interference effects are a consequence 
of the temporal phase coherence of the wavefunction. The ability to vary this coherence by 
varying the temporal separation of the avoided crossings renders twisted rapid passge with 
adjustable twist strength into a temporal interferometer through which qubit transitions 
can be greatly enhanced or suppressed. Cubic and quartic twist were explicitly considered 
in this paper as they are, respectively, the simplest examples of odd-order and even-order 
polynomial twist in which these interference effects are expected to occur. Although we have 
focused on these two cases, we do not mean to suggest that these pulses are the best of all 
possible twisted rapid passage pulses. A search is currently underway for other twisted rapid 
passage pulses that might produce stronger destructive interference, and hence, faster, more 
fault tolerant quantum CNOT and NOT gates (see below for further discussion). We have 
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FIG. 13: The transition probability P{t) for non-adiabatic rapid passage with quartic twist with 
A = 5.0 and 774 = -4.6 x lO""^. 

seen that this interference mechanism can be used to pump qubit energy-levels, as well as 
to strongly quench qubit transitions during non-adiabatic twisted rapid passage. Although 
cubic twist proved to be more effective at pumping than quartic twist, quartic twist was 
found to be much more effective at quenching qubit transitions. We have seen that quartic 
twist allows qubit inversion to be done both non-adiabatically and at fidelities that exceed 
the threshold for fault tolerant operation. The marriage of operational speed with reliability 
is a direct consequence of the destructive interference that is possible between the 3 avoided 
crossings that can arise during rapid passage with quartic twist. 

B. Implementing Quantum CNOT Gate 

We now describe a procedure for implementing a quantum CNOT gate using twisted 
rapid passage in the context of liquid state NMR. If the liquid has low viscosity, one can 
ignore dipolar coupling between the qubits, and if the remaining Heisenberg interaction 
between the qubits is weak compared to the individual qubit Zeeman energies, it can be 
well-approximated by an Ising interaction Under these conditions, the Hamiltonian (in 
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FIG. 14: The transition probability P{t) for non-adiabatic rapid passage with quartic twist with 
A = 5.0 and 774 = 1.6 x 10^^. Note the shghtly reduced vertical scale. 



frequency units) for the control (c) and target (t) qubits is: 



(24) 



Here uj^ {ujt) is the resonance frequency of the isolated control (target) qubit; J is the Ising 
coupling constant; and Uc > uJt > t^J ■ We choose the single qubit computational basis 
states (CBS) to be |0) = | t) and |1) = | i). Thus the two-qubit CBS are: |00) = | TT); 
|01) = I ti); |10) = I []); and |11) = | ||), and they are the eigenstates of Hd- The 
energy-levels (in frequency units) are shown in Figure El where 



(25) 



Given this energy-level structure, we can implement a quantum CNOT operation on the 
iwo qubits by sweeping through the resonance using twisted rapid passage. Decoupling 
18^ is used to switch off the dynamics of the control qubit so that only the target qubit 
responds to the rapid passage pulse. Since the two states |00) and |01) are not resonant, 
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FIG. 15: The transition probability -P(t) for adiabatic rapid passage with quartic twist with 
A = 0.5 and r/4 = 6.45 x 10"^ 
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FIG. 16: Energy-level structure appropriate for implementing a quantum CNOT operation using 
twisted rapid passage. The corresponding energies (in frequency units) appear to the right of the 
energy levels. 



they do not respond to the twisted rapid passage pulse. Thus, 

|00) ^ |00) 
|01) - |01) . 



(26) 



TABLE II: Transition probabilities for quartic twist with A = 5.00 and r/4 in the range [3.95 x 10 ^, 
4.04 X 10-^]. 



774 ( X 10 j 




3.95 


2.0 X 10^^ 


3.96 


1.3 X 10-2 


3.97 


6.8 X 10-3 


3.98 


3.6 X 10-3 


3.99 


9 X 10-^ 


4.00 


4 X 10-^ 


4.01 


8 X 10-^ 


4.02 


3.9 X 10-3 


4.03 


1.0 X 10-2 


4.04 


1.7 X 10-2 



On the other hand, for the |10) and |11) states, the combination of decoupHng and sweeping 
through the a;+ resonance means that only the target qubit has its spin flipped. Thus, 

|10)^|11) , , 

(27) 

|11)-|10) , 

and we see that this procedure implements a quantum CNOT operation on the two qubits. 



C. Experimental Realization 

Because of the fundamental sig nificance of quantum CNOT gates to quantum computing 



and quantum error correction |l4 



15 



, 3| 1 it is hoped that the feasibility of using rapid pas- 
sage with quartic twist to implement this gate might be tested experimentally (see penulti- 
mate paragraph of Section|l}. Experimental realization of polynomial twist (/)n(t) = {2jn\BV^ 
should be possible through an adaptation of the procedure used by Zwanziger et. al. jsf to 
realize quadratic twist. Thus: (1) the driving rf-field is linearly polarized along the x-axis in 
the lab- frame with ^^.(t) = 2b cos (f)rf{t)', (2) the resonance offset at (see eqn. (^Q)) is produced 



by linearly sweeping the detector frequency ujdet{t) through the resonance at the Larmor fre- 
quency ojq such that Udetit) = luq + {2at/Ti)] and (3) twist is introduced by sweeping the 
rf-frequency Urfif) = (prf through the resonance at uq in such a way that uj^fif) = u^et — <Pn- 
It is worth noting that the resonance condition uJrf{t) = ujq is identical to our existence 
condition for avoided crossings, eqn. 0. Note that in our paper the external field inversion 
takes place over the time-interval (—To/2, To/2); the external field crosses the x-y plane at 
t = and is initially aligned along the — z direction. The Appendix provides a translation 
key which relates the theoretical parameters of this paper to the experimental parameters 
of the Zwanziger experiments 0,0]. 

Before leaving the subject of experimental realization of rapid passage with quartic twist, 
two further remarks are in order. First, to insure that all qubits are inverted when a spread 
of resonance frequencies occurs, it is necessary to require that the frequency sweep cover a 
large enough interval that the entire spread of resonance frequencies is included in it. This 
gaurantees that all qubits will have passed through resonance by the end of the frequency 
sweep. Second, a range of rf field strengths can also be accomodated so long as aT/2 ^ bmax- 
This condition insures that the frequency sweep begins far from resonance for all rf field 
strengths, and that transitions will continue to occur only near the avoided crossings. One 
therefore anticipates that in this case also, the interference effects will continue to occur as 
predicted. For reasonably good samples, magnets, and rf sources these constraints can be 

satisfied, and the interference effects presented above should be readily observable. This is 

1 

in fact what is found experimentally j^. 



D. Other Pulses 



Having introduced twisted rapid passage with polynomial twist, and pointed out the 
possible advantages of quartic twist for quantum computing, it is natural to ask how quartic 
twist compares with the more familiar vr-pulse which can also be used to implement quantum 
CNOT and NOT gates. We begin by comparing the inversion time for quartic twist with 
that of a comparable vr-pulse. We focus on quartic twist with A = 5.00 and 774 = 4.00 x 10~^ 
as this choice of parameters yields a gate fidelity F = 0.99996 (see Section IIVC|) which 
exceeds the threshold for fault tolerant gate operation Fft ~ 0.9999. We now show that this 
case achieves fault tolerant operation while simultaneously matching the inversion speed of 



a vr-pulse. In the notation of Ref. the basic experimental parameters for twisted rapid 
passage are A, B, ui, and T, and they are related to our theoretical parameters by eqns. flA8|) 
and ()A3|) . T continues to denote the duration of the twisted rapid passage pulse. Because 
a twisted rapid passage sweep must begin far from the avoided crossings, A and ui cannot 
be chosen independently. In the rf-frame, the asymptotic effective magnetic field must lie 
near the z-axis so that tan 6* = ooi/A ~ 0.1. Choosing uJi = 4000 Hz gives A = 4 x 10^ Hz. 
Both of these values can be achieved with existing NMR technology. Writing / = Ui/A, 
eqn. ()A9|) gives 



With A = 5.00, this gives 

T4 = 2 msec . 

By comparison, the inversion time for a 7r-pulse with rf-amplitude Ui = 4000 Hz is = 
Tc/ui = 0.8 msec. Thus, twisted rapid passage with quartic twist is clearly capable of 
matching the inversion speed of a comparable vr-pulse while still exceeding the threshold for 
fault tolerant operation. On the other hand, the error probability for a typical vr-pulse is 
P ~ 10^'^ due to, for example, inhomogenities in the rf field amplitude. This corresponds 
to a fidelity F ~ 0.999 so that, unlike the equally fast quartic twist pulse which acts fault 
tolerantly, the vr-pulse falls short of the threshold for fault tolerant operation Fth ~ 0.9999. 

We hope in the future to examine higher order versions of polynomial twist to determine 
whether they have more effective quenching and/or robustness properties than cubic and 
quartic twist. We have also done preliminary work on the interesting case of periodic twist: 
(f)(t) = TTpsinut. As we have seen, polynomial twist only allows 1-3 avoided crossings to 
occur during rapid passage. One can show that periodic twist allows the number of avoided 
crossings that occur during rapid passage to be modified through variation of the twist 
amplitude p and frequency u. We intend to explore how the interference effects considered 
here are modified when more than 3 avoided crossings can occur. 
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APPENDIX A: CONNECTION BETWEEN THEORY AND EXPERIMENT 



For ease of comparison with Refs. \^ and 0, we choose F{t) = — 6cos0„(t)x — 
6sin0„(t)y + atz in eqn. (^. The Hamiltonian in the detector frame is then 

H(t) at b ^ . , / N 

—— = — (T^ - - cos (pn{t) cr^- T sm (p„{t) ay . 
n n n h 

Here 0„(t) = {2/n)Bt"', and to avoid confusion with the notation of Ref. P], we have 

switched the symbol used for the twist strength in the main body of this paper: B ^ B. 

Transformation to the rf-frame is done using the unitary operator U{t) = exp[— (i/2) 0„(t) cr^] 

so that H{t) H{t): 



H{t) 



n 



b 

h 

2b 



¥ T 

2at 

and I = a/2. 

The experimental Hamiltonian in the rf-frame appears in eqn. (12) of Ref. j^: 



(Al) 



h 



Comparing eqns. ()A1|) and ()A2|) gives 



2b 



(A2) 



(A3) 



and 



Integrating eqn. ()A4|) gives 



2at 



<l>rf{t") 



dt" 



T/2 



h 



T 



T 



(A4) 



(A5) 



In this paper, we have parameterized time such that t'" G [ —T/2, T/2], and T is the duration 

of the twisted rapid passage pulse. Defining 

_ t'" 1 
r- Y + 2 ' 

it follows that r G [0, 1]. Introducing t' = t'" + T/2 and t = t" + T/2, eqn. (jMI becomes 



<Prf{t) 



dt' 



2aT 



2i3T" 
T 



(A6) 



rf = 4>det - 4>n] with (t)det = ^^0 + 



As explained in the caption of Figure 2 of Ref. P], 
2A{t — 1/2); and generalizing to polynomial twist, 0„ = nB^r — l/2)"~^/T, where B is the 
symbol used in Ref, 

into 4>rf = 4>det ^ 



for the twist strength. Plugging these expressions for 0det and 
6„, and integrating gives 



dt' 



Equating eqns. ()A6|) and ()A7|) gives 



2A 



nB 



T 



n-1 



aT 



B 



2B 



n 



rpn 



Using eqns. ()A3|) and ()A8|) in the definition of A (see discussion following eqns. 

4U| 



A 



Using eqns. ()A3|) . ()A8|) and eqn. (fT5|) with n = 3 and 4 gives 



^3 = 7 



3 5cJi 



and 



4 
Ecu? 



2A^T 



(A7) 

(A8a) 
(A8b) 
gives 
(A9) 

(AlO) 
(All) 



respectively. The results of this Appendix give the connection between our theoretical 
)arameters and the experimental parameters A, B, ui, and T of the Zwanziger experiments 
^, 0]. In Section 13 these formulas are used to calculate the inversion time for a twisted 

rapid passage pulse with quartic twist. 
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